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Abstract 



We discuss the line shape for an internal state transition for bosonic atoms 
confined as a one-dimensional gas. Typical line shape is an edge singularity 
due to the absence of Bose-Einstein Condensation in such systems. 
PACS numbers: 03.75.Fi 

Recent progress in trapped bosonic gases has stimulated a lot of activities in the study 
of many-body effects in such systems. Most theoretical papers concentrate on effects due 
to the existence of the Bose-Einstein Condensate (BEC) and the associated macroscopic 
wave-function. With current techniques, it is likely that one can as well study an effectively 
one-dimensional (ID) trapped gas This situation occurs when the bosons are put in 
a trap which is tightly confining in two directions ( y and z) but is essentially free in the 
third (x), and with particles occupying only the lowest quantized subband for motions in 
the y-z plane. This system is unique in that, for any finite interaction among the bosons, 
there is no Bose-Einstein Condensation even at zero temperature. p[. The long- wavelength 
quantum fluctuations destroy the phase coherence of the system in the macroscopic limit. 
For a given density n, the occupation number Nq of the lowest single particle state, though 
in general increases with the size and hence the total number of particles in the system, is 
nevertheless only a neglig ible fraction of the total {^^^^^ ). Though many properties 
of ID bose gas have been studied in the past 0, the subject remains somewhat academic 
since there exists no physical realization. The possibility of realizing this ID bose gas has 
already stimulated some more theoretical work on this system recently 

The present ID bosonic system is closely related to its fermionic counterpart; and they are 
collectively known as the "Luttinger liquids" 0. One-dimensional fermionic systems have 
been studied much more extensively in the context of quantum wires, carbon nanotubes 
etc. There the one-dimensional nature leads to the absence of fermi liquid behavior, so 
the properties of the system is qualitatively different from its non-interacting counterpart. 

The atomically trapped ID bose gas offers us an excellent new opportunity to study the 
peculiar properties of quantum systems in reduced dimensions, in particular the consequence 
due to the absence of BEC. Not only that this is the first ID bosonic system available, it 
provides advantages over the mentioned ID fermionic cases in that the present system is 
intrinsically clean, and moreover numerous atomic and optical tools can now be used. It 
may even be possible to further study novel systems such as mixtures etc. 

In this paper we consider one such example, namely the absorption line shape of an 
optical transition in this ID bose gas. We imagine initially a quasi- Id system of bosons in 
their internal ground state (referred to as the 'a-atoms' below). An incident optical wave 
excites one of the atoms to a different internal state. We are interested in the absorption 
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line shape of this process; i.e., the probabihty that the absorption takes place at a given 
frequency Q. We shall consider the problem at T = 0, paying particular attention to the 
special feature due to the absence of BEC. An analogous experiment has already been done 
in a three dimensional trap for the Is 2s transition in hydrogen [^. For a uniform 3D bose 
gas at T = 0, the line is expected to be narrow, Lorentzian (see, e.g. 0, c.f. however [Q) 
with a width governed by the small 'gaseous parameter (s)' (^3^0^)^/^. Here n^d is the particle 
number density in 3d and a('s) is (are) the s-wave scattering length(s) for scattering among 
the bosons. For this three dimensional case, the main weight of the transition comes from 
exciting a particle from momentum p = 0. There is a macroscopic number of such particles. 
For a weakly interacting gas, the potential felt by the rest of the bosons is only slightly 
modified after the transition. We shall see that the situation in ID is very different due 
to the absence of BEC. Exciting a single particle necessitates a substantial rearrangement 
of the relative motion among the rest of the particles, i.e., emission of a large number of 
phonons. The basic line shape is typically an edge singularity. 

We shall then consider an atomic trap with tight confinement in the y-z directions. For 
definiteness we shall assume that the confinement potentials in these directions are isotropic 
and harmonic, with frequency uj±. We shall consider the case where the atoms are limited 
within < X < L but otherwise free. A weak harmonic trap potential along the x-axis has 



been considered in the past [|IO],|TT|, and the low dimensional effects have been shown to be 
significantly reduced. We shall assume that such a potential is absent here. We shall further 
ignore the effect due to finiteness of L, a condition which we shall return below. 

The bosonic system is described by a Hamiltonian containing the kinetic energy and a 
short-range interaction among the bosons. Anticipating that we shall eventually study the 
system in the quasi-ld limit, the field operator for the original bosons, referred to hereafter 
as the a-bosons, is expanded as V'a(^) = i'a{x)xQ{r±) + J2j 4'aj{x)xj{r^) where ipa{x) is 
the annihilation operator for the a-bosons in the lowest subband (0) where the transverse 
wavefunction is Xo(ri)- i^aji^Ys ( j = 1,2,...) are the corresponding operators for the 
higher subbands (wavefunctions Xj{^±))- the ground state \G >, ipaj{x)\G >= since all 
particles reside in the lowest subband. We shall see that ipaj for j = 1, 2, ... will not appear 
anywhere below, and can pretend that the expansion of ipa{'^ will consist thus only the term 
involving the lowest subband. The extent of the wavefunction xo iii the y-z direction will 
be denoted as a±. For harmonic trapping potentials in these directions = 



We shall not write down the interaction part of the Hamiltonian involving the a-bosons 
explicitly since we won't need it below. For delta-function type interactions it is possible to 
write down all the formal many-body wavefunctions |T^. Unfortunately manipulations of 



such wavefunctions are typically very mathematically involved. Moreover we shall eventually 
be interested in 'final' states where an a-atom has been excited internally. Such wavefunc- 
tions, where effectively there is a 'foreign' atom present, are not known (with an exception 
noted below). We thus proceed rather differently and follow Haldane 0, concentrating only 
on the low energy excitations which are density oscillations in the system. The field operator 
i^ai^) ( = ^pO-pe^^^ I \fV) describing the motion along x is re-expressed in terms of density 
n{x) and phase 0(x) by i^ai^) = [n(x)Y^'^e'^'^^^\ The effective Hamiltonian for the Id motion 
can be written as 



Ho = ^ldx 



K 



(1) 
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where is related to the number density fluctuations Sn{x) by = nSn. and 
9{x) must be considered as operators satisfying the commutation relation Vx'0{x')] = 

i'K5{x — x'). V is the density (sound) velocity wave of the system, and K is the (dimen- 
sionless) Luttinger liquid interaction parameter. vK = imoh/m and K is related to the 
compressibility of the system via K = Trvh{dno/dfi) where Uo is the linear number density 
and fi is the chemical potential of the system. K depends on the interaction among the 
bosons. For short range interactions, K in principle have already been obtained 0JT^ but 
there is no general analytic form known. We simply note here that K = 1 for the strongly 
repulsive (impenetrable) limit while ^ cxd if the interaction among the bosons is weak. 

Hq can be diagonalized easily |^ since it is quadratic. It is useful to introduce bosonic 
operators bq and 6j which describe the sound modes of the system: 

'/>(^) = -^Ei^r/v«^(&^-M (2) 

Then, apart from a constant, Hq = J2q^o ^^qb\bq where the mode frequencies are uiq = v\q\. 

Now we consider the external optical fleld responsible for the transition in inter- 
nal state. Since the atom with a new internal state is distinguishable from the origi- 
nal atoms, we shall call the resulting atom the 'c-atom'. For deflniteness we consider 
only the 'Doppler-free' part of the spectrum ( no momentum transfer from the exciting 
laser beam(s)) and assume that the perturbation responsible for the excitation is uniform 
throughout the ID bose gas. The relevant part of the Hamiltonian can be then written as 
Hex = w J d^r {ip\{r)'il)a{r) + '0l('^V'c(?^) where ipd'f^ is the fleld operator for the atom in 
the excited internal state. We are interested in the rate of transition as a function of the 
excitation frequency Vt. This rate is given simply by the golden rule: 

liSl) = 2tiY.\< F\H,.,\G > \^5in - (Ep - Eg)) (3) 

F 

where \F > are the set of flnal states. To proceed further we need to know the fate of the 
c-atom. There are many possibilities and we shall simply consider two extreme limits: 

(1) An atom in the internal state c is not affected at all by the trapping potential for 
the atoms in state a. This is possible if, e.g., this potential is due to a laser with frequency 
near a dipole resonance for the a- atoms but far away from any of those for c. In this case 
the c-atom escapes from the trap and no longer interacts with the remaining a-atoms. 

(2) The c-atom also feels the strong transverse conflnement potential. It remains inside 
the trap and continues to interact with the a-atoms. 

We shall treat these two cases in turn. 

Case 1: In this case the flnal state is simply a product state, \F >= |$ > x\p>, with 
the — 1 a-atoms in state |$ > (with x momentum —px) and the c-atom in the plane wave 
state \p > with momentum p. The rate for transition into states with given p is 

Ij,{n) = 27r\fpf I < > - {E^ + f/2m + u- Eg)) (4) 
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where ff=J (P'r±e~^^^'^^xo{'^±) is a form factor, and u is the difference in internal energy 
between the a- and c- atoms. This rate is thus proportional to the spectral function for 
annihilation of a particle in the ID bose gas (hereafter Px —>■ P ) 

B{p,uj') = J2\< '^K\G > - {E^ - Eg,n-i)) (5) 

with Lu' = Q + Eg — Eq^n-i — p^/2m — u and where Eq^n-i is the ground state energy of 
— 1 atoms in the ID trap. B{p,uj') is given by —^lmGi{p,uj') where 

Gi{x, t) = i< V-KO, 0)V',(x, t) > hit) (6) 

The required Green's function can be evaluated by finding the long distance and time 
behavior of Gi{x,t). Substituting eq (@) into eq(|) one finds Gi ~ (a;2^^2f2-)i/4y h{t). 
The Fourier transform Gi{p,uj') and hence B{p,uj') for a given p has the scaling form 
~ h{uj' — f |j9|)/(c(j'^ — v'^p'^y^"' with a = Since 1 < K < oo the line shape is of the 
form of an edge singularity. For the three dimensional bose gas, the corresponding B{p,uj') 
is proportional to a delta function at the Bogoliubov mode frequency (with an additional 
small incoherent part). The lack of a delta function in the present case reflects the absence 
of a condensate. If p can be independently measured this experiment would be analogous 
to angle resolved photoemission spectroscopy employed extensively in recent studies of high 
T-c cuprates [0. 

Case 2: In this case the final result depends on the interaction between the a and c 
atoms, which has the general form gad dripl{^'4'c(yf^4'i{f^4'a{r) ■ Here gac is defined in 3D 
and is related to the 3D scattering length a^c by gac = 4:nhaac/Tn . We express ipcir) also in 
subbands ipcir) = J2j i^cj{x)xcj{rj_). The final states in eq (|^) in general can involve terms 
that are off-diagonal in the band index j for the c-atom due to the interaction between a 
and c. Such processes correspond to the possibility that the transverse motion of the c- 
atom be modified due to the interaction with the a-atoms. For a tight trap however, these 
contributions are small in the dilute limit. They involve the parameter nogac/'^(^\^L ~ 
aac/l- Here / = l/rio is the average interparticle spacing among the a-atoms. Ignoring thus 
these contributions, we have = Y^j \fj\'^Ij{Q), i.e., the transition line then becomes a 
superposition of 'lines' involving transition to final states where the c-atom is within a given 
j-th subband, with a weight given by the form factor fj = J df±[xo{'i^±)xtj{'^±)]- The shape 
of each of these lines is determined by Im[Dj^(n)] where 

Df(t) = -iJ dxi J dx2 < i^iixi,t)ipcj{xi,t)i;lj{x2)tpa{x2) > h{t) (7) 

which we shall now compute. The effective Hamiltonian is given by the sum of Hq involving 
only the a-atoms and Hcj which describes the motion of the c-atom within its j-th subband 
and its interaction with the a-atoms. Hq has already been given in eq(0) and the effective 
Hcj is given by H^^ + if™* where 

consists of the 'one-body' (kinetic, trap and internal energy) terms and 
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H'c' = 9ac,n I dx4Ax)i^,,{xWa{x)Ux) (9) 



is due to the interaction between the a and c atoms. Here is the subband energy for the 
c-atom in its j-th subband and Qacjj = dad (P'i^XcjXo{f^±)\'^ depends on j. For small j it is 
of order Qac/ T^a\ if the trap potentials for a and c atoms are similar. 
It is convenient to express Dj^(t) in the 'interaction' picture: 

Df(t) = -ij dx,j dx2 < ^l(xi,t)7A,,(xi,t)fe-*/o^cf(0^*' ^1^.(^2)^^(3,2) > /,(i) (10) 

where T is the time-ordering operator and the expectation value is to be calculated with the 
Hamiltonion Hq+H^^^ . One can rewrite eq(|TDp in momentum space and expand the exponent 
involving the interaction, if™* then produces 'scattering' terms where the momentum of the 
c-atom is changed from one value to the other. The resulting calculation does not seem to 
be tractable analytically in general. The problem however can be simplified significantly 
if we assume that gac,jj is small (condition given below) and concentrate again on the low 
frequency limit, i.e., just above the threshold for transition into the j-th subband (given by 
^cj = ^cj — ^0 — fi' + u + gacjjnid to lowest order in gac)- In this case transitions must be made 
to states where the momenta p of the c-atom are small. Provided that p << mv, then the 
kinetic energy (of motion along x) transferred to the c-atom p^ /2m is much less than that 
to the density oscillations of a-atoms vp. Thus so long as we are interested in frequency 
deviations Af2 from the threshold which satisfy Af2 << mu^, one can ignore the kinetic 
energy of the motion of the c-atoms along x. For an impenetrable bose gas this amounts to 
limiting ourselves to Af2 << Tc'^h'^ /ml'^ . Ignoring thus the first term in H^j^ and rewriting 
the result back in real space, we obtain 

Bf{t) = -iLe-'^^^' < V^l(0,t) re-'/o^'".-*(*')'^*>„(0) > h{t) (11) 

where the effective interaction Hamiltonian Ha^int acts only on the a-atoms and is given by 
fl'ac,jjV'a(0)V'a(0). This rcsult can be understood physically as follows: the external optical 
field produces an a — >■ cj transition at t = and at a general location say Xq, annihilating 
an a-atom while producing a c-atom there in the j-th subband. Since the c atom moves 
with velocity p/m whereas the sound waves move with the speed f , at small p the c-atom 
is essentially stationary. For t > 0, due to the c-atom created, a delta-function interaction 
at xq acts on the a-atoms. We are interested in the overlap between the initial state with 
an a-atom destroyed at Xq and the final states with this extra interaction potential. This 
overlap is independent of the location of the transition xq which then has been set to 0. The 
problem has thus become similar to that of X-ray absorption in solid state physics. There 
initially ( t < 0) one has an electron gas in its ground state. At t = an X-ray photon is 
absorbed which creates a charged nuclei with an extra electron added to the electron gas, 
which at t > feels an extra local potential due to the charged nuclei. Notice however there 
are slight differences. Here for t < we have an equilibrium ID interacting bose gas and at 
t > there is one less boson than initially. 

The correlation function in (pA]) and thus the line shape can be found analogous to the 
X-ray problem [Q. The expectation value needed can be rewritten as < ?/'''' (0)e~*^*?/'(0) > 
where H = Hq + Hint- Hq = J2 ft-^qblbg and Hint = gac,jjSn{0) can both be written in 
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terms of bq and The annihilation operator ip ~ e**^ acts like a displacement operator 
for these bosons: e'"^bqe~'^^ = bq — 1 2gLA' 1"^^^- Thus the expectation value is the same as 
< Tq'"^' So^9ac,jj-T^V]^)Sn{o,t )dt ^ calculated with the Hamiltonian Hq. Here we have defined 
vn = v/K. The calculation can be easily done since Hq is quadratic in the boson operators bq 
and The long time dependence is given by, apart from a part oscillating sinusoidally with 
t which contributes to the line shift, Df{t) ~ and thus the line shape ~ l/[Afi]^~"j 

with a'j = jxl^ ~ m^]^" Notice that the line shape is in general different for transition to 
different subbands of the 'c-atom'. 

It is instructive to compare this result with that of the X-ray problem when the interac- 
tion between the nucleus after the X-ray and the electrons is modelled by a delta-function 
interaction of strength V. There the long time behavior of the corresponding correlation 
function is given by and the line shape is ~ l/[Af2]-^~"'^ with ax = [1 + N{0)V]'^, 

where iV(0) is the density of states at the fermi level for the electron gas fl^. Here, the 
interaction strength is replaced by gacjj and the density of states is replaced by |^ = 
In both cases the power law line shapes (for ax 7^ and a'^ 7^ 0) are due to 'orthogonal 
catastrophe', that the overlap between the state just after the excitation and the new ground 
state of the system vanishes. There are however two differences: the -|- sign in ax is replaced 
by a — sign since now there is one less boson at t > rather than one more electron in the 
X-ray absorption problem. The decay of D^{t) in time is slower {a'j is smaller) if Qac > 0: 
near the location of the excitation, the reduction in local density of the a-atoms is also what 
an interaction with Qac > prefers. The factor 1/2K arises from the fact that we have a 
one-dimensional interacting bose gas rather than a three-dimensional fermi liquid. 

As the repulsive interaction gacjj between the a and c atoms increases, a'j decreases 
and the line sharpens up. The above quantitive result for the exponents a'j however holds 
only for small Qacjj ( << TTVNh ). In the limit where interactions among all the bosons are 
strongly repulsive, one reaches the impenetrable limit where the bosons cannot pass each 
other. In this limit the statistics, i.e., the indistinguishability among the a-bosons and the 
distinguishability between the a- and c- bosons become irrelevant. One can check that the 
wavefunctions written down in [|r^ for a system of identical bosons only can be generalized 
to the present case if we replace one of the bosons by the foreign c-atom. It can be further 
verified that the absorption line shape becomes delta functions (one for each subband cj) in 
this impenetrable limit. 

The experimental observation of the line shape discussed above should be feasible. The 
possibility of obtaining a quasi-lD bose gas has already been discussed in ref [|I|. The above 
line-shape is applicable so long as Af2^^ = ui. For hydrogen atoms and with / ~ 1/xm, 
Ui ~ lOOkHz, much larger than the present available experimental resolutions. It is sufficient 
for the temperature T to be small compared with the line- width. The condition T « uji 
is satisfied for T < ^K, a temperature readily achievable. A finite length L of the system 
will break the line into a set of discrete sub-lines with separations of order vr^^/mL^, but 
the shape of the 'line' can still be observable as long as L » I. 

I thank Ite Yu for a useful correspondence. 
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